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Abstract 

We study the behavior of the interquark potential in lattice gauge theories at high temper- 
ature, but still in the confining phase, and propose a new observable which could play in this 
regime the same role played by the Liischer term in the low temperature limit. This quantity 
is related to the exponent of the power prefactor in the effective string partition function or, 
equivalently, to the coefficient of the logarithmic correction in the interquark potential and, 
as for the usual Liischer term, its value does not depend on the particular gauge group under 
consideration or on the form of the effective string action used to model the flux tube. In this 
respect it can be considered as a universal signature of the effective string behavior of the flux 
tube. As a test of our proposal we studied this quantity with a set of high-precision numerical 
simulations in the (2+1) dimensional SU(2), SU(3) and SU(4) Yang-Mills theories and in the 
Z2 gauge model, always finding a perfect agreement with the predicted values. 



PACS numbers: 12.38.Gc, 11.15.Ha, ll.10.Wx, 11.15.Pg, 12.38.Aw 



1 Introduction and motivation 



It is by now rather well understood that the flux tube joining a static quark- antiquark pair in 
the confining regime of a generic lattice gauge theory (LGT) can be described by an effective 
string theory pQ. The most relevant consequence of this description is the presence of a term 
in the interquark potential proportional to 1/r (where r is the interquark distance) with a fixed 
coefficient (the well known "Liischer term" pQ) which can be predicted from the theory and is 
not affected by corrections due to higher order terms in the effective string action. In this sense 
the Liischer term can be considered as a universal signature of the effective string behavior of the 
flux tube and in the past years it was the first goal of any numerical study on the subject. More 
recently it has been realized that this term is only the first of a series of universal corrections in 
the perturbative expansion of the effective string action in powers of or 1 [2j[3] opening the way 
to more refined numerical test of the effective string picture (for an updated review of lattice 
results, see [1]). 

While the first numerical tests, more than twenty years ago, were performed using Wilson 
loops in these last years it has become more frequent to choose a finite temperature setting 
and look at the interquark potential using Polyakov loop correlators (or torelons, depending to 
the choice of the space and time directions on the lattice). These finite temperature tests can 
be performed in two very different regimes. One of them is the "open string channel", which 
in the LGT language is the low temperature regime: it corresponds to lattices with a very long 
size in the "inverse temperature" direction and hence very long Polyakov loops. In this regime 
the dominant effective string correction to the potential is again the Liischer term and higher 
order (possibly universal) string corrections appear as contributions with higher powers in 1/r 
in the interquark potential. The other possible choice is the "closed string channel", which in 
the LGT language corresponds to the high temperature (but still in the confining region) regime. 
In this case the size of the lattice in the compactified direction is just above the deconfinement 
length and the Polyakov loops are much shorter than the typical interquark distance. In this 
regime the dominant effective string correction is linear in r and has the effect of decreasing the 
string tension which becomes temperature dependent. As the temperature increases the lattice 
size in the compactified direction becomes shorter and shorter and accordingly the string tension 
decreases and at the end it vanishes leading to the deconfinement transition. From a physical 
point of view, this high T regime is particularly interesting since in this limit the behavior of the 
flux tube and of its quantum fluctuations could strongly influence the transition from hadrons 
to deconfined quarks as T c is approached. The fact that the Liischer term gets mixed with 
higher order corrections makes it difficult to use it as a universal signature of the presence of an 
effective string. Indeed in this regime one usually has to compare the numerical results with the 
expectation of the whole effective string action which is not universal anymor^]. 

Also in this regime, it would be interesting to have a universal signature of the effective string 
which could play the same role played by the Liischer term in the low temperature ( "open string" ) 
limit. As we shall show in detail below, the natural candidate for such a role is the coefficient 

1 Indeed in this limit the typical question is how accurate are the results from the Nambu-Goto model, which is 
the simplest possible effective string theory and in several cases turns out to fit the numerical data very well [4] . 
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of the logarithmic correction in the interquark potential (or, equivalently, the exponent of the 
powerlike prefactor in front of the effective string partition function). This term, which grows 
logarithmically with r, is not affected by higher order contributions of the effective string action 
and only depends on the nature of the conformal field theory (CFT) to which the effective string 
model flows in the infrared limit. In particular, the presence of this term is a signature of the free 
bosonic nature of the quantum fluctuations of the flux tube in this limit i.e. the fact that they 
can be described by the c = 1 CFT of a free uncompactified bosonic degree of freedom. In this 
respect it can be considered as a high temperature analogue of the Liischer term and, exactly as 
the Liischer term, it is universal and has a simple linear dependence on the number of transverse 
dimensions. 

The main goal of the present paper is to study this observable and test its universality in a 
set of LGTs both with continuous and discrete gauge groups. In this respect it is important to 
note that a projection of the Polyakov loops onto their zero transverse-momentum component 
(which is normally performed in this type of lattice computations [3]) would wash away such a 
term. Therefore in our computation we did not perform a projection onto the zero transverse- 
momentum component, and the correlators we investigated also involve higher momenta. This 
demands a special computational effort, and a combination of state-of-the-art simulation and 
error-reduction algorithms, with large statistics. 

As we shall see, in all the cases we studied (which are certainly described by different effective 
string actions, since their critical behavior in the vicinity of the deconfinement point is different), 
the coefficient of this logarithmic term in the potential is always the same and it agrees with the 
predicted value. This makes this observable a robust signature of the effective string behavior in 
the closed string channel and a good candidate to play the role of the Liischer term in the high 
T regime of confining LGTs. 

This paper is organized as follows. In section [2j we remind the construction and the most 
important properties of the bosonic effective string model, discussing, in particular, the Nambu- 
Goto action. Section [3] recalls the basic features of Yang-Mills theories in 2 + 1 dimensions and 
introduces the setup of our lattice simulations, whose results are reported in section 0J Section [5] 
presents a discussion, and some concluding remarks. 

2 Effective string theory 

In the confined phase of Yang-Mills theories, the zero-temperature ground state interquark po- 
tential Vq(t) for a pair of fundamental, static color sources QQ at asymptotically large distances 
r is linearly rising with r: this can be interpreted as evidence for the formation of a stable "tube" 
of chromoelectric flux lines between the two charges, with an energy proportional to its length. 
In the infrared limit, such tube can be interpreted as a string-like object, with negligible intrinsic 
width and no transverse structure. If there is one such string in the confining Yang-Mills vacuum, 
then the translational and rotational symmetries of spacetime get spontaneously broker]! down to 
the subgroup of translations parallel to and rotations around the string world sheet. This leads 

2 The breaking is actually explicit, if the sources are not located at infinity. 
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to the existence of D — 2 massless Nambu-Goldstone bosons (the transverse fluctuation modes of 
the string) , which are expected to be the relevant degrees of freedom of an effective theory, valid 
at energy scales well below the intrinsic mass gap of the confining theory. 

It is important to emphasize that this effective string picture is expected to provide a universal 
low-energy description, i.e., one independent of the underlying gauge theory (as long as it is 
confining) and of its microscopic dynamics details. 

For an open string of finite length r, the quantization of the transverse vibration modes 
(treated as non- interacting) leads to a Casimir effect, which manifests itself as a 1/r terrrjl in the 
ground-state interquark potential [I]: 

Vo(r) ~ a r - ^=^, (1) 

to a tower of evenly spaced excited levels: 

E n (r) ~ V (r) + -n, n G N, (2) 
r 

and to a logarithmic broadening of the mean square width w 2 of the fluctuations in the midpoint 
of the string, as a function of the string length: 

D — 2 

w2 ( r ) ~ t; ln ( r / r o) (3) 

Zttij o 

(where tq denotes a fixed length scale). 

In order to describe the complete partition function for the sector of the gauge theory con- 
taining a pair of static color sources, one needs to generalize these predictions by including string 
interactions, which are encoded in the complete effective string action S eS . During its time-like 
evolution in a Euclidean setup — possibly at a non-negligible temperature T — , each static source 
traces out a straight worldline which is a Polyakov loop winding around the compact time di- 
rection. The two-point Polyakov loop correlation function G(r, T) is then written as the string 
partition function: 

G{r,T) = (V* (r)V(6)) = J Vhe~ s ^ h \ (4) 

where the functional integral of the transverse fluctuations of the string h (which are the degrees 
of freedom in the physical gauge) is to be performed over connected world-sheet configurations 
having the Polyakov loops as their boundary. 

The functional form of S cB is not known a priori, but it can be expanded as a series of 
derivatives of h, including both bulk and boundary terms. At each given order in the number 
of derivatives and fields, the various terms can be constrained by Poincare symmetry and open- 
closed string duality [3] (see also [E]). Somewhat unexpectedly, it turns out that the first few 
orders in this series exactly reproduce the expansion of the Nambu-Goto action |13j . which is just 
proportional to the area of the string world-sheet: 

Sng =<tq [ d 2 A. (5) 



Note that such a term has nothing to do with a short-distance Coulomb term. 
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The effective string expansion shows deviations from eq. ((5|) only at high orders [3 [16], which 
are at the limit or beyond the resolution of current lattice calculations. In particular, for the 
open string channel in D = 2 + 1 dimensions, the leading-order corrections (which are induced 
by boundary terms appearing in the effective action) are predicted to be proportional to r~ 4 , 
and have been recently investigated in high-precision simulations [IB]. For the closed string 
channel, the leading-order corrections appear at an even higher order, and are proportional to 
r -7 . This justifies considering the Nambu-Goto action as a very accurate approximation of the 
actual effective string action, within the precision of present numerical simulation results. 

Thanks to its particularly simple form, the energy levels and properties of the Nambu-Goto 
action can be easily derived analytically; in particular, the open string spectrum reads |17j : 




E° n = a rJl + - „-___ , n € N. (6) 



It is well-known that the Nambu-Goto string is generically affected by an anomaly, leading to 
breakdown of rotational invariance (except for the critical number of spacetime dimensions D = 
26 [IB]), and that, being non-polynomial, it is non-renormalizable. However, these issues are 
not of concern to us here, given that in the present context one just regards eq. (0) as (an 
approximation to) the action for an effective theory valid at low energies only and that the effects 
of the anomaly in non-critical dimensions are suppressed in the long-string limit [19] . 

Approximating S eS by S NG on the r.h.s. of eq. the Polyakov loop correlation function in 
D = 2 + 1 dimensions at a finite temperature T can be expressed as a series of modified Bessel 
functions. This result was obtained a few years ago by Liischer and Weisz [2] using a duality 
transformation and then derived in the covariant formalism in [14] . In d = 2 + 1 dimensions one 
finds a tower of Kq Bessel functions: 

oo 

G(r,T) = J2cnK (rE° n ) (7) 

n=0 

with arguments involving the Nambu-Goto energy levels for a closed string: 



K= U ^h^(n-hl new (8) 



T y cr V 24 

In particular, in the large-distance regime rT S> 1, eq. ([7]) is dominated by the term depending 
on the lowest energy level: 

G(r,T)~K (rE%). (9) 
Following ref. [20J, one can also define a temperature-dependent string tension a = <r(T) as: 



a = E ^T = aJl-^-, (10) 



which is predicted to vanish at a critical temperature T C NG = ^2>ao/i:. 
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Given that the two-point Polyakov loop correlation function G(r, T) can be related to the 
finite-temperature heavy-quark potential V(r, T) via: 



V(r,T) 



T\nG(r, T) 



(11) 



and that the modified Bessel function K n {z) has the following expansion: 



K n (z) 




—z 



1 + 



4n 2 - 1 16n 4 - 40n 2 + 9 



8z 128z 2 



+ 0(z- 3 ) 



(12) 



for large values of its argument z, one gets: 



V(r, T) ~ rTE c + - In (rE%) + 



T 



T 



(13) 




where we dropped an irrelevant additive constant, and neglected terms which are suppressed by 
higher powers of (rT) . 

At this point, a remark is in order: the corrections to the effective string action can be written 
in a series of powers of (err 2 ) -1 , where a is the string tension at the temperature T, and vanishes 
in the T — > T~ limit, for a second-order phase transition. This imposes a further constraint on the 
validity range of our expansion, and, in particular, implies that the model cannot be expected to 
hold all the way up to T c . In our analysis, we shall include only data corresponding to err 2 > 1.5, 
which — based on previous studies — appears to be a reasonable criterion. 

Within the r and T ranges satisfying these conditions for the validity of the expansion, eq. (fT3|) 
predicts thus a characteristic, logarithmic contribution to the confining heavy-quark potential. 
It is important to stress that, for the values of r and T usually studied in lattice simulations, 
this logarithmic term is much larger than the statistical errors. Moreover, as mentioned in 
the introduction, it is universal in the sense that any effective string model which has a free 
bosonic CFT as its large-r limit must contain this contribution which can thus be considered as a 
distinguishing feature of this class of effective string models. At the same time it is important to 
stress that such a term would be washed away, if one — as it is often done in lattice computations — 
were to project the Polyakov lines to their zero transverse momentum components. From the 
point of view of the numerical lattice computation, this makes the observation of the effect more 
challenging, and requires the use of error-reduction techniques which are discussed in the next 
section. In the remaining part of this paper we shall mainly focus on this contribution. 

3 Confining gauge theories in 2 + 1 dimensions and their lattice 
regular izat ion 

In this section, we first introduce SU(iV) Yang-Mills theories in 2 + 1 dimensions in subsection 13. 11 
then we discuss their lattice regularization in subsection 13.21 as well as the Z2 lattice gauge 
theory (subsection 13. 3ft . Finally, we define the observables evaluated in our lattice simulations in 
subsection 13.41 
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3.1 Yang-Mills theories in 2 + 1 dimensions 



SU(iV) gauge theories in D = 2 + 1 spacetime dimensions are particularly interesting models 
characterized by non-trivial dynamics and sharing many qualitative features with ordinary Yang- 
Mills theories in D = 3 + 1. They can be defined in terms of the following Euclidean functional 
integral: 

Z = J VAe~ s \ 5 E = | d^TrF^, (14) 

where the bare square gauge coupling g$ has the dimensions of a mass. Thus the natural expansion 
parameter for bare perturbation theory calculations at a momentum scale k is the dimensionless 
ratio g$/k, which leads to a non-trivial infrared structure |21j . Similarly to Yang-Mills theories 
in D = 3 + 1 dimensions, these theories are asymptotically free at high energy; one difference 
w.r.t. the four-dimensional case, however, is that Yang-Mills theories in D = 2 + 1 dimensions are 
superrenormalizable (rather than renormalizable) , i.e. they have a finite number of ultraviolet- 
divergent Feynman diagrams. 

At low energies, these theories are linearly confining with a finite mass gap and a discrete 
spectrum. At zero (and low) temperature, the physical states are color-singlet glueball states 
which can be classified according to the irreducible representations of the 0(2) group and charge 
conjugation properties (for A > 2). 

Like in the D = 3 + 1 case, SU(A) Yang-Mills models in D = 2 + 1 dimensions undergo a 
deconfinement transition at a finite critical temperature T c , at which the Zjy center symmetry of 
the low-temperature vacuum gets spontaneously broken by a non-vanishing expectation value of 
the Polyakov loop in the thermodynamic limit. Again, similarly to what happens in D = 3+1 [23j . 
the deconfinement transition is of second order for "small" gauge groups, and turns into a more 
and more strongly first-order one when the gauge group size is increased — in agreement with the 
intuitive picture of a more and more "violent" transition when a larger number of gluons get 
liberated [24J. In particular, the transition in D = 2 + 1 dimensions is of second order for the 
gauge groups SU(2) and SU(3), while it is a weakly first-order one for SU(4), and a stronger 
first-order one for SU(A > 5) [2~5lJ26| (whereas in D = 3 + 1 the transition is of second order 
only for the SU(2) gauge group). Here we focus our attention on the cases of N = 2, 3 and 4 
colors, and study these theories by regularizing them on a lattice, as discussed in the following 
subsection. 



3.2 Lattice regularization 



We focus our attention on the cases of N = 2, 3 and 4 colors and study these theories non- 
perturbatively by regularizing them on a finite, isotropic cubic lattice A of spacing a and volume 
V = Li x Lt = (Ng x Nt)a 3 . The lattice dynamics is defined by the standard Wilson lattice 
gauge action: 



E 

xgAl<a</3<3 



1 



1 

iv J 



ReTr U n 



with: (3 



2N 

~ 2 ' 



(15) 
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and: 



U af3 (x) 



U a {x)Up(x + aa)Ul(x + a/3)ul(x). 



(16) 



Periodic boundary conditions are imposed all directions, so that the shortest lattice size Lt is 
related to the temperature of the system via: T = 1/Lf. The functional integral appearing in 
eq. (|14p is then replaced by the finite-dimensional multiple integral: 



3 

2l= HI] dU a (x)e 

•J A 1 



cE 
L 



(17) 



xGA a=l 



(where dU a (x) denotes the Haar measure for the U a (x) link matrix), while expectation values of 
gauge-invariant physical observables O are obtained from: 



<e>> 



3 

y- 1 n n du ^) ° 



qE 
L 



(18) 



and are estimated numerically by Monte Carlo sampling over a set of {U a (x)} configurations. In 
the following, we denote the number of configurations as n conf . In our simulations, the samples of 
thermalized, independent configurations were produced using a 1+3 combination of heat-bath |27| 
and overrelaxation updates [28] on SU(2) subgroups [29j , 

Although high-order lattice perturbation theory computations are also available, typically 
numerical lattice simulations (both in D = 3 + 1 and in D = 2 + 1 dimensions) are run in a 
regime of strongly coupled dynamics, where the determination of the spacing a as a function of 
the bare gauge coupling (or, equivalently, of /3) is done non-perturbatively, extracting a through 
the lattice determination of some reference quantity relevant for low-energy scales (such as, for 
example, the zero-temperature string tension o"o, or the critical deconfinement temperature). 
The systematic uncertainty related to the choice of a physical observable to set the scale has 
quantitatively modest impact 0(a 2 ). In our computations, we used the scale determination from 
ref. [6] for SU(2), from ref. [15] for SU(3), and from ref. [26] for SU(4), yielding: 



<7oa ~ < 



1.324(12) . 1.20(11) 

3 + ^^~ 

3.37(1) , 3.90(25) , 50.1(1.8) " 
P p 2 1 p 3 

6.491(68) , 3.7(3.2) , 533(37) 
P JP 1 



for SU(2), 
for SU(3), 
for SU(4). 



(19) 



Just to get an intuitive picture in terms of scales relevant for real-world QCD, eq. (|19p can be 
translated into a scale for a in fm, by defining gq to be equal to (440 MeV) 2 . 
Table Q] summarizes the basic technical information about our simulations. 



3.3 Z 2 lattice gauge theory in three dimensions 

As mentioned in the Introduction, we also compared the results from Yang-Mills theory with those 
from the Abelian Z 2 gauge theory in three dimensions: this is a lattice model whose fundamental 
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gauge group 


N'f x N t 


P 


a a 2 


a 


physical volume 


T/T c 


^conf 


SU(2) 


no 2 x 8 


9.0 


0.0262(1) 


0.072 fm 


8.7 2 x 0.58 fm 3 


3/4 


5.6 x 10 4 


SU(3) 


120 2 x 8 


21.34 


0.0293(1) 


0.077 fm 


9.2 2 x 0.61 fm 3 


3/4 


8 x 10 4 


SU(4) 


120 2 x 8 


40.10 


0.0282(6) 


0.075 fm 


9.0 2 x 0.60 fm 3 


3/4 


4 x 10 4 



Table 1: Parameters of our SU(iV) lattice simulations (see the text for the definition of the various 
quantities). 

degrees of freedom are cr a (x) variables defined on the oriented bonds of a three-dimensional cubic 
grid with periodic boundary conditions and taking values in (the fundamental representation of) 
the cyclic group of order two: Z2 ~ ({1,-1},-). In analogy to eq. (|15|). the dynamics of this 
model is defined by the action: 

ffgaugc = -/3ga U gc &ap{x), with: cr a/3 (x) = a a (x) 073 (x + aa ) a a (x + a/3) 073 (x) , (20) 

xGA l<a<,9<3 

which is invariant under Z2 gauge transformations flipping the signs of the six a a (x) variables on 
the bonds touching the site x. 

Despite its deceivingly simple definition, this model is characterized by highly non-trivial 
dynamics and a rich phase structure. For values of /3 gaugo less than /3 gaugej crit . = 0.76141346(6) [30] 
(and for sufficiently large lattice sizes), the model is in a confining phase, while above that 
critical value, it is in a deconfining phase. The bulk phase transition separating the two phases is 
a second order one, which enables one to investigate an effective continuum limit of this confining 
lattice model by approaching /3 gaugc , crit. from below (by virtue of the fact that ooa 2 — > for 
/5 g au g o — > /3~ ugCi crit .) an d scaling the lattice size appropriately. Furthermore, in the confining 
phase there also exists an infinite-order roughening transition at /3 gaugo , rough . = 0.47542(1) [31], 
separating the /3 gaugc < /3 gaugo , r OUg h. regime where strong-coupling expansions hold, from a "rough" 
phase /3 gaU g C , rough. < /3gau g o < /^gaugo, crit. in which the theory can support massless excitations. 
Furthermore, when this model is defined on a lattice with a compact dimension of size Lj, one 
can observe a finite-temperature deconfinement transition when Lj becomes shorter than the 
inverse of the critical temperature. 

An interesting feature of the Z2 gauge model defined by eq. ([20]) is that, through a duality 
transformation [32], it can be mapped exactly to the tridimensional Ising spin model with s(x) € 
Z2 variables, which is defined by: 

Sspin = -fispfa ^ s(x)s(y), with: sinh(2/3 apin ) sinh(2/3 gaugo ) = 1, (21) 

where (x, y) denotes that the sites x and y are nearest-neighbors. This duality implies a one-to- 
one mapping between the free energy densities of the two models in the thermodynamic limit, 
and allows one to express expectation values of gauge-invariant observables in the Z2 gauge 
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model as ratios of partition functions of the spin system, with a set of appropriately chosen 
antiferromagnetic couplings, leading to a major efficiency boost in simulations of this model. 

In particular, in this work we studied the %2 g au g e model analyzing two ensembles of configu- 
rations corresponding to two different temperatures T = T c /2 and T = 2T c /3, taken from ref. [8]; 
technical details are summarized in table [2J 



Ensemble I 


Ensemble II 


iV| x N t 


96 2 x 12 


iV 2 x N t 


96 2 x 9 


/^gauge 


0.746035 


/^gauge 


0.746035 


cr a 2 


0.018943 


do a 2 


0.018943 


a 


0.062 fm 


a 


0.062 fm 


physical volume 


5.92 2 x 0.74 fm 3 


physical volume 


5.92 2 x 0.55 fm 3 


T/T c 


1/2 


T/T c 


2/3 


^conf 


3.2 x 10 9 


^conf 


3.2 x 10 9 



Table 2: Parameters of the Z2 simulations analyzed in this work. 



3.4 Observables 

In order to investigate the logarithmic term appearing in eq. (|13p . we estimate numerically the 
expectation values of (on-axis) two-point correlation functions of Polyakov loops at a finite tem- 
perature T = l/(aN t ) < T c : 

G h (r,T) = y [UU dU a (x)V*(ar,0)V(0,0) e~ s * , (22) 

with: 

Nt-l 

P(an x , an y ) = Tr [ U t {an x ,an y ,an t ). (23) 

n t =0 

To measure G L (r, t), we used a highly efficient multi- level algorithm [TU], which allows one to 
achieve exponential error-reduction. The error analysis was done using the standard jackknife 
method [33], taking cross-correlations into account. 

Following |11] and [6] it is particularly convenient to introduce the following quantities: 

^' T > - l^l&T^TY |24) 

A{rT) _ g ^ + ^-.T) (25) 

Note that, in the continuum limit a — > 0, Q(r,T) tends to the first derivative of V(r,T) w.r.t. r: 

dV 

limQ = — , (26) 
a->0 or 
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so that it can be interpreted as a lattice version of (minus) the interquark force. On the other 
hand, A(r, T) is a dimensionless quantity proportional to the discretized derivative of the force: 

r 2 d 2 V 

limi = --— . (27) 
a ^o T dr 2 y ' 

It is easy to see that these quantities are the finite temperature version of the observables 

introduced in [11]. In particular Q(r, T) coincides in the low-T limit with the "force" F(R) of [11] 

while A(r, T) is related to the "central charge" c(R) of [llj as follows 

A(r,T) = ^c(r). (28) 

Using eq. ([9]) , eq. (PTOl) and eq. (fT2j) we may estimate the large-r limit of these two observables 
in the framework of the Nambu-Goto effective string model: 

In these two functions the universal logarithmic term in which we are interested is encoded in 
the coefficient \/{2r) of the term proportional to T in eq. (j29[) and in the constant 1/2 in eq. (|30p . 
This last value, thanks to the normalization of A(r,T) in eq. (I25p . exactly coincides with the 
coefficient of the logarithmic term in the interquark potential. In the next section we shall use 
our Monte Carlo simulations to estimate these quantities in the SU(iV) and Z2 models that we 
studied. 



4 Results 

4.1 General remarks on the fitting procedure 

Let us briefly describe a few features of our fitting procedure. 

• As previously discussed, the effective string description holds for large enough interquark 
distances. Following previous works on the finite temperature interquark potential, we 
decided to include in the fit only values of the interquark distance r > r min such that 
a"(T)r 2 in ~ 1.5. To fulfill this constraint we had to choose r min = 9a for SU(iV) models, 
r min = 10a for Ising, Ensemble I and r min = 12a for Ising, Ensemble II. As a consistency 
check, it is interesting to notice that with this choice in all the cases the reduced x 2 turned 
out to be always of order unity. 

• In the SU(iV) case, the data for Q turned out to be strongly cross-correlated and thus we 
had to perform our fits taking the cross-correlation matrix into account. It is important 
to stress that a naive fit to the data neglecting cross-correlations would lead not only to 
different values of the reduced x 2 an d of the statistical uncertainties but also to differences in 
the best fit values larger than the quoted errors. On the contrary, the Ising data for different 
values of r are completely uncorrelated since they were obtained in different simulations. 
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4.2 Analysis of the SU(2) data 



Table [3] reports our numerical results for the "discretized" interquark force Q(r,T) denned by 
eq. (|24p . as obtained from SU(2) simulations at fixed T = 3T c /4. For convenience, we express 
both r and Q in units of (the appropriate power of) the lattice spacing a. 



r/a a 2 Q 


r/a a 2 Q 


r/a a 2 Q 


2 0.037433(46) 

3 0.030958(56) 

4 0.027600(64) 

5 0.025553(72) 

6 0.024154(84) 

7 0.023118(94) 


8 0.02232(11) 

9 0.02170(12) 

10 0.02117(12) 

11 0.02072(13) 

12 0.02034(15) 

13 0.02000(15) 


14 0.01971(16) 

15 0.01949(18) 

16 0.01926(19) 

17 0.01906(20) 

18 0.01892(22) 

19 0.01876(24) 



Table 3: Results for a 2 Q(r,T), as a function of the interquark distance r/a, from the SU(2) 
simulations at T = 3T c /4. 

We fitted the values of Q(r,T) according to 



a 2 Q(r,T)\ T=3Tc/4 = s + — + , 
' r/a (r/a) z 

and found the following best fit values for the parameters 



(31) 



s = 0.01530(37) b = 0.0668(58) c = -0.087(27) 

with a reduced Xr = 0.75. 

The universal correction in which we are interested is encoded in the parameter b which 
according to the analysis discussed in the previous sections should be given by 

aT 1 

»= T = „=<»*». 

which turns out to be in remarkable agreement with our result. 

This is further confirmed by the analysis of the A(r, T) values (which can be easily obtained 
from the data reported in table EJ). We fitted these values to 



m 



finding 



A(r,T)\ T=3Tc/i = k — 

' r/a 



k = 0.528(28), m = -1.09(28), with xL = 1.6, 



(32) 



which is again in perfect agreement with the expected value k = 1/2. 

From the first fit we can extract the value <t(T) = 0.01530(37) for the finite temperature string 
tension at aT = 1/8. Using the value a^a 2 reported in eq. (|19fl . we may obtain a "Nambu-Goto" 
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prediction for a(T) using eq. (fl~0[) . which turns out to be a 2 a NG (aT = 1/8) = 0.01605(6), at two 
standard deviations from the observed value. This indicates, as already observed in [22], that for 
the (2+1) SU(2) LGT the Nambu-Goto string represents a rather good approximation but, with 
the new generation of high precision data, small deviations start to be detectable. 

Finally, in the same way, we may obtain predictions for the subleading corrections in the 
two fits. We find for the c term in the first fit c NG ~ —0.1216(5) and for m in the second 
fit m NG — —1.946(8). Both values are of the same order of magnitude and of the same sign 
of the results, but they are not compatible within the errors. This discrepancy agrees in sign 
and magnitude with the analogous deviations from the Nambu-Goto ansatz observed in [22] and 
summarized in the coefficient C3 evaluated in |22| . 

A similar pattern is observed in the SU(3) and SU(4) cases, which we describe in detail below. 

4.3 Analysis of the SU(3) data 

For the SU(3) gauge group, our results for a 2 Q(r, T), as a function of r/a, are reported in tabled 



r/a a 2 Q 


r/a a' 2 Q 


r/a d 2 Q 


2 0.042353(20) 

3 0.035013(24) 

4 0.031241(28) 

5 0.028971(33) 

6 0.027447(37) 

7 0.026351(42) 


8 0.025525(47) 

9 0.024857(54) 

10 0.024326(61) 

11 0.023876(68) 

12 0.023501(76) 

13 0.023184(88) 


14 0.022881(96) 

15 0.02262(11) 

16 0.02238(12) 

17 0.02223(14) 

18 0.02196(15) 

19 0.02178(17) 



Table 4: Results for a 2 Q(r,T), as a function of the interquark distance r/a, from the SU(3) 
simulations; the temperature is fixed at T = 3T c /4. 



In this case, the fit to eq. (|29j) yields 

a 2 Q(r,T) | T=3Tc/4 = 0.01884(44) + 



0.0612(74) 0.063(34) 



, with xL = 1.6, (33) 



r/a ( r / a ) 

to be compared with 0.01946(6), 0.0625 and —0.1003(5) from the Nambu-Goto model, while the 
fit of A(r, T) gives 



A(r,T)| T=3Tc/4 = 0.529(46) - L °° ( * 5 } , with X L = 0.9, 



(34) 



to be compared with k = 1/2 and m NG = —1.605(6). 
4.4 Analysis of the SU(4) data 

We report our results for Q(r, T) for SU(4) gauge group in table [5j they can be fitted to eq. ([29 

0.0634(70) 0.063(32) 



as: 



a Q(r, T)| T=3Tc /4 = 0.01721(43) + 



r/a 
12 



(r/a) 



, with X ? ed = 1.4, (35) 



to be compared with 0.0183(4), 0.0625 and —0.107(2) from the Nambu-Goto model, while the fit 
of A(r, T) gives: 

A(r,T)\ T=3Tc/4 = 0.470(30)-^^, with X L = 0-7, (36) 

to be compared with k = 1/2 and m NG = -1.71(4). 

It is interesting to observe that, although the disagreement with the Nambu-Goto values 
increases as N increases, the expectation value of k (or equivalently of b in the first fit) is in 
perfect agreement with the expected value k = 1/2 (or b = 1/16), thus showing that this term 
is indeed universal (in the same sense in which the Liischer term is universal) and is a robust 
indicator of the free bosonic nature of the transverse excitations of the effective string. 



r/a a 2 Q 


r/a a 2 Q 


r/a a 2 Q 


2 0.040834(29) 

3 0.033492(33) 

4 0.029726(37) 

5 0.027468(42) 

6 0.025953(47) 

7 0.024871(53) 


8 0.024044(58) 

9 0.023379(65) 

10 0.022847(73) 

11 0.022403(81) 

12 0.022030(91) 

13 0.02170(10) 


14 0.02142(11) 

15 0.02118(13) 

16 0.02093(14) 

17 0.02068(16) 

18 0.02052(18) 

19 0.02037(20) 



Table 5: Results for a 2 Q(r,T), as a function of r/a, from our SU(4) simulations at T = 3T c /4. 
The results for all the three gauge groups, and the fitted curves, are shown in fig. (TJ 



4.5 Analysis of the Z 2 data 

As usual, the data in the Ising case show a slightly different behavior. 

In the ensemble I (which has T/T c = 1/2) the agreement with Nambu-Goto is remarkably 
good, while for the ensemble II (which is closer to the critical point, with T/T c = 2/3) sizeable 
deviations from the Nambu-Goto predictions can be observed. On the contrary, the values of 
k and b are fully compatible, like for SU(iV), with the effective string expectations for both 
ensembles. 

Our results for Q(r,T) at T = T c /2 are reported in tabled and can be fitted to: 

■Aj(r, r)| TrfV2 _ 0.01485(2) + ^> - , with xL = 0.7, (37) 

to be compared to the prediction from the Nambu-Goto model: 0.01487, 0.04167 and —0.58 
(due to the precision in the determination of o"o, errors in the Nambu-Goto estimates of these 
parameters can be neglected), while for A(r,T) the corresponding result^! read: 

A(r, T)\ T=Tc/2 = 0.47(3) - ^j^-, with X L = 0.85, (38) 

4 Note that, since in the Z2 simulations we determined Q(r,T) only for even values of r/a, in this case we 
modified the definition of A(r, T) in eq. ()25[) by replacing a with 2a. 
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Figure 1: Left-hand-side panels: Results for the "discretized interquark force" a 2 Q(r,T) defined 
in eq. (|24p . as a function of the source separation in units of the lattice spacing r/a, at fixed 
T = 3T c /4, together with the fitted curves, for SU(2) (top), SU(3) (middle) and SU(4) (bottom). 
Right-hand-side panels: Same, but for A(r, T) as a function of r/a, together with the fitted curves. 
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to be compared with k = 1/2 and m = —1.40. 



r/a a 2 Q 


r/a a 2 Q 


r/a a 2 Q 


8 0.0192898(84) 
10 0.0185110(94) 
12 0.0179654(93) 
14 0.017568(11) 
16 0.017242(11) 


18 0.017010(13) 
20 0.016811(13) 
22 0.016633(13) 
24 0.016500(13) 
26 0.016361(14) 


28 0.016285(15) 
30 0.016179(16) 
32 0.016125(15) 
36 0.015960(17) 
40 0.015838(18) 



Table 6: Results for a 2 Q{r,T), function of the interquark distance r/a, from the % 2 simula- 
tions at T = T c /2. 

Similarly, table [7] displays the results for Q(r,T) for the ensemble at T/T c = 2/3; in this case 
the results of the fits are: 

, s 0.0522(40) 0.076(34) , 9 
a 2 Q(r,T)\ T=2Tc/3 = 0.01137(11) + - V with X L = 1.8, (39) 

for Q(r, T) (for which the coefficients predicted by the string model are 0.01067, 0.0556 and 
—0.145, respectively), and: 

A(r,T)\ T=2Tc/3 = 0.47(6) - ^P, with Xr 2 ed = 1.1, (40) 

for A(r, T), to be compared with k = 1/2 and m = —2.60. 



r/a a 2 Q 


r/a a 2 Q 


r/a a 2 Q 


8 0.016872(14) 
10 0.015904(15) 
12 0.015194(17) 
14 0.014728(18) 
16 0.014320(19) 


18 0.014021(21) 
20 0.013755(21) 
22 0.013614(21) 
24 0.013460(23) 
26 0.013292(25) 


28 0.013131(25) 
30 0.013022(25) 
32 0.012888(26) 



Table 7: Results for a 2 Q(r,T), as a function of the interquark distance r/a, from computations 
for the Z2 model at T = 2T c /3. 

Figure [2] displays our results from the two Z2 ensembles at T = T c /2 and T = 2T c /3, together 
with the fitted curves. 

It is interesting to observe that (as in the SU(iV) case) the values of k and b are fully compatible 
with the effective string expectations for both ensembles, even though in the second ensemble the 
other parameters disagree within the errors from the Nambu-Goto estimates. 

Finally, let us mention that, as already noticed in [HE], the best fit values obtained in this 
second ensemble turn out to agree with a truncation of the Nambu-Goto model up to the 0(T 4 /cjo) 
term in the expansion of eq. (|10p . which gives s = 0.01138. 
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Figure 2: Same as in figure [TJ but for the Z2 gauge model at T = T c /2 (top) and at T = 2T c /3 
(bottom). 
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5 Discussion and conclusions 



The results presented in the previous section show that, in all cases that we studied, the numerical 
data are in good (but not exact) agreement with the predictions from the Nambu-Goto string 
model, eq. (|29p and eq. (|30p . This is particularly true for the SU(iV) models, while in the Ising 
case one observes increasing deviations as the temperature increases. 

However in all the cases, even when the Nambu-Goto prediction for the string tension is not 
compatible within the error with the first term of eq. (|29p . the coefficients of the logarithmic 
term in the potential (which is encoded in the b term of eq. (J29]) or, equivalently, in the k term 
of eq. (|30p ) always agree with the predicted value. In particular, the dimensionless ratio A(r,T) 
at fixed T clearly tends to the expected value k = 1/2 in the limit of large quark-antiquark 
separations r — > oo. 

The fact that A tends to the asymptotic value 1/2 indicates that, whatever the exact effective 
string action describing the world-sheet dynamics in this temperature range, its large-r limit is 
represented by a c = 1 conformal field theory describing a free uncompactified bosonic degree of 
freedom. In this respect, this logarithmic term can be considered as a high-temperature analogue 
of the Liischer term and, exactly as the Liischer term, it is universal (as we have seen in our 
simulations) and has a simple linear dependence on the number of transverse dimensions. 

As mentioned in the introduction, the Liischer term in the high temperature regime, due to the 
peculiar modular symmetry of the string partition function, becomes a correction proportional 
to the interquark distance. The same happens to all the higher order corrections due to the 
higher order terms in the effective string action. All these terms sum up to give the temperature 
dependence of the string tension and as a consequence the Liischer term cannot be resolved 
anymore. On the contrary, all higher-order terms of the effective string action do not contribute 
to the logarithmic correction, which thus maintains its free bosonic value k = 1/2, irrespective of 
the details of the effective string action, and thus it can be used as a universal signature of the 
effective string description. 

As mentioned above, the precision of our data is sufficient to reveal clear quantitative devi- 
ations from the predictions of the Nambu-Goto string. These are most clearly exhibited by the 
coefficients of the 0(r~ 2 ) terms in the fitted expressions for Q, eqs. (|3"T|) . (j33j) . (|35|) . (|37l) and 
(j39"P . which agree in sign with the expectations from eq. (f29j) . but are different in amplitude. Our 
data analysis shows that such deviations are not compatible with statistical fluctuations: rather, 
they reveal a real physical effect, namely that in this temperature range the idealization of the 
confining flux tube as a Nambu-Goto string already starts to fail. 

This is not surprising and is consistent with similar observations already pointed out in the 
literature (for an updated review see, for instance, ref. [1]). In particular, ref. [7] also found that, 
for SU(3) Yang-Mills theory in D = 2 + 1 dimensions, sizeable deviations from the Nambu-Goto 
model show up for temperatures higher than 0.7T C — although the "temperature-dependent string 
tension" remains close to its Nambu-Goto expectation in that regime. Our results are in full 
agreement with these findings, and generalize them to groups other than SU(3). 

Concerning the gauge group dependence, one interesting observation is that our results for 
the Z2 gauge model show somewhat more pronounced deviations w.r.t. the Nambu-Goto model. 



17 



In particular, the constant term in eq. (|39p tends to deviate from the corresponding predictions 
for the temperature-dependent string tension a 2 a. This agrees with previous observations [EEE], 
that a truncated Nambu-Goto prediction is actually (accidentally) in better agreement with the 
numerical data than the full expectation. On the other hand, the analogous results for the SU(2), 
SU(3) and SU(4) gauge groups are in better agreement with the Nambu-Goto prediction, and do 
not reveal a particular dependence on the number of colors. 

On the theoretical level, another observation is that, if one squeezed the lattice size in the 
direction orthogonal to the plane defined by the two Polyakov loops (i.e. the direction in which 
the effective string fluctuates), then at some point, when the lattice size becomes of the same 
order of the effective flux tube width, the string would start to wind around the compactified 
transverse direction (we assume periodic boundary conditions in all the lattice directions). At 
this point it would be more natural to model the flux tube fluctuations in terms of a compactified 
bosonic string theory. In this case, the leading term in the large-distance behavior of A would 
vanish; this effect could be used for precision numerical tests of the flux tube width. 

To summarize, in this work we have presented the results of a high-precision numerical inves- 
tigation of confining gauge theories in D = 2 + 1 dimensions, at finite temperature. We studied 
the behavior of the force between static sources in the fundamental representation, and compared 
the results of our lattice simulations for SU(2), SU(3) and SU(4) Yang- Mills theories, as well as 
for the Z2 gauge model, with the predictions of an effective theory based on a confining bosonic 
string described by the Nambu-Goto action. In particular, we focused on an observable which can 
be interpreted as a finite-temperature analogue of the Luscher term, and which can be obtained to 
high precision from correlators of Polyakov lines — provided one does not perform a projection to 
their zero transverse- momentum component. This poses a technical difficulty, which we overcame 
by using advanced simulation algorithms, including multi- level error-reduction methods |10| . 

We found good agreement between the string theory model and lattice results, and a very 
mild dependence on the gauge group rank in SU(iV) Yang-Mills theories. This confirms the 
strong similarities between these theories, which has already been remarked in other works, both 
in D = 2 + 1 and in D = 3 + 1. We also found clear quantitative deviations from the Nambu- 
Goto predictions, which show up in the subleading correction terms in the interquark force and 
its derivative. Remarkably enough, even in the cases in which we observe deviations from the 
Nambu-Goto expectation, the coefficient of the logarithm term in the potential keeps the expected 
value k = 1/2 and can be used as a precise signature of the effective string in a finite temperature 
geometry. 
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